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Abstract Comparisons of solar magnetic-field measurements made in different 
spectral lines are very important, especially in those lines in which observa- 
tions have a long history or(and) specific diagnostic significance. The spec- 
tral lines Fei 523.3 nm and Fei 525.0 nm belong to this class. Therefore, this 
study is devoted to a comprehensive analysis using new high-precision Stokes- 
meter full-disk observations. The disk-averaged magnetic-field strength ratio 
R = 5(523. 3)/B(525.0) equals 1.97 ± 0.02. The center-to-limb variation (CLV) 
is i? = f .74 — 2.43/i+3.43/i^, where fi is the cosine of the center-to-limb angle. For 
the disk center, we find R = 2.74, and for near-limb areas with fi = 0.3, R equals 
f .32. There is only a small dependence of R on the spatial resolution. Our results 
are rather close to those published three decades ago, but differ significantly 
from recent magnctographic observations. An application of our results to the 
important SOHO/MDI magnetic data calibration issue is discussed. We conclude 
that the revision of the SOHO/MDI data, based only on the comparison of 
magnetic field measurements in the line pair Fei 523.3 nm and Fei 525.0 nm 
(increasing by a factor of 1.7 or 1.6 on average according to recent publications) 
is not obvious and new investigations are urgently needed. 
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1. Introduction 

Many tasks of solar and solar-terrestrial physics urgently require precise, quanti- 
tative information about the distribution of magnetic fields across the solar disk. 
However, the extremely complicated spatial structure of magnetohydrodynamic 
parameters in the solar atmosphere leads to differences in magnetic field data 
from different spectral lines. This is the reason why comprehensive analyses of 
magnetic-field measurements made in different spectral lines are important for 
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diagnostics of solar magnetism and for the analysis of different data sets. In the 
first case, it is better to use observations performed at the same instrument, to 
avoid possible influences of different instrumental, method, time-dependent, and 
other effects. In the second case, the aim is to construct composite long-term 
data sets combining different observatories to study, for example, the temporal 
variations of solar magnetic parameters on different time scales. 

A very significant part of solar magnetic-field observations has been performed 
(and are performed up to now) in the spectral line Fei 525.0 nm. Therefore, 
comparisons of observations in other spectral lines with observations in this 
particular line are extremely important. The atomic parameters of Fei 525.0 nm 
(large Lande factor, g = 3, small lower level excitation potential, xi = 0.12 eV) 
stimulated vigorous discussions about the reliability of measurements in this 
line, which are not finished yet. Different authors (some aspects of this issue 
are given in the discussion in Dcmidov et ai, 2008 and in the recent paper of 
Ulrich et ai, 2009) offer different correction factors to reduce raw measurements 
in Fei 525. Onm to the "true" ones. 

Because of its low temperature sensitivity, one of the lines used as a "stan- 
dard" for such correction is the rather strong line Fci 523.3 nm with joff = 1-3 
and Xi = 2.93 eV (Harvey and Livingston, 1969; Howard and Stcnflo, 1972; 
Frazier and Stenflo, 1972; Frazier and Stenflo, 1978; Uhich, 1992; Uhich et ai, 
2002; Tran et ai, 2005; Ulrich et ai, 2009). However, the combination of this 
line and Fei 525.0 nm is far from being perfect for diagnostics of solar magnetic 
fields, because they have different thermal dependencies, and Fei 523.3 nm is 
not intrinsically narrow (see Solanki, 1993). With much more success, it could 
be used as an example for the thermal line ratio. As shown by Socas-Navarro 
et al. (2008), reliable information about the magnetic-field parameters in the 
quiet regions (which cover most of the solar surface) can be achieved only using 
a combination of several lines with particular characteristics. 

Nevertheless, observations in the discussed lines arc still important in the 
context of the cross-calibration of different data sets, especially for the data 
recalibration problem of the Michelson Doppler Imager (MDI) onboard the Solar 
and Helio spheric Observatory (SOHO), see http://soi.stanford.edu/magnetic and 
Solar News (issues No. 19, 2007 and No. 25, 2008). The SOHO/MDI observations 
of magnetic fields are performed in the spectral line Nil 676.8 nm with poff = 
1.43. The main reasons for two subsequent serious revisions of the magnetic field 
data from MDI, which are widely used in the scientific community, are based 
on the comparison of observations in Fci 523.3 nm and Fei 525.0 nm made at 
the Mount Wilson Observatory (MWO) and published by Tran et al. (2005) and 
Ulrich et al. (2009). According to the first paper, based on the results of Ulrich 
(1992) and Wang and Shcclcy (1995), observations in Fci 525. Onm must be 
multiplied by the factor = B(523.3)/B(525.0) = 4.5-2.5 sin^ 6 or, according 
to the second paper, by the factor 6^^ — 4.15 — 2.82 sin^ 9, where 9 is the ccnter- 
to-limb angle. As a consequence, to fit SOHO/MDI data to MWO observations 
in Fei 525.0 nm corrected in this way, they must be multiplied by a factor that 
is a function of 9, see Figure 6 and Table 3 of Tran et al. (2005) and whose mean 
value over the disk is 1.7 . According to Ulrich et al. (2009), the factor is 1/0.619 
= 1.615 for 9 < 30°. 
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The fact, that the new MWO results strongly contradict the previous ones of 
Howard and Stenflo (1972), who found that the ccnter-to-limb variation (CLV) 
of the coefficient is S"-^ = 0.48 + 1.33/i, with /i — cos 6*, requires independent 
measurements. For disk center, the factor is only 1.8 instead of 4.12 according 
to Uhich et al. (2009). Their factor between B(523.3) and 5(525.0) can be up 
to 5.5 (see their Figure 3), which can hardly be explained by "the fact that 
the shifted Zeeman components at A525.0nm are shifted beyond the sampling 
passband of the MWO Babcock magnetograph" (Ulrich et al., 2009). 

The above mentioned circumstances make an additional exploration of solar 
magnetic field observations in Fei 523.3 nm and Fei 525.0 nm an urgent task. 
The entire spectral-line profile should be considered. This is the main aim of our 
study. 



2. Observations and Analysis 

All previous observations concerning the relationship between magnetic-field 
data in Fei 523.3 nm and Fei 525.0 nm were obtained with photoelectrical mag- 
netographs. Drawbacks of such data are evident, because direct information 
about the distribution across the spectral line profile of the Stokes parameters 
(Stokes V/Ic in the case of longitudinal field) is missing - the only reliable indi- 
cator of magnetic field in the spatial range of line formation. Some assumptions 
and calibrations arc necessary to obtain magnetic-field strengths from intensity 
variations in the exit slits of the magnetograph, which might be different for 
different spectral lines in the case of simultaneous observations. 

Spectropolarimetric observations, as used in the present study, are much 
more reliable and informative. To study the problem of the comparison of solar 
magnetic-field measurements in combination of these lines, several series of ob- 
servations in different regimes have been performed with the Stokes-meter of the 
Solar Telescope for Operative Predictions (STOP) at the Sayan Solar Observa- 
tory (SSO). Basic information about this instrument and methodical issues can 
be found in Demidov et al. (2002) and Demidov et al. (2008). STOP is equipped 
with a linear, 29 mm wide Toshiba TCD CCD detector with 3424 pixels (height 
200 ^m, width 8 /^m), which allows us to obtain high-precision measurements 
of the Stokes parameters / and V/Ic in several spectral lines simultaneously. 
The solar surface is scanned following a 2D-raster, and the spatial resolution 
depends on defocussing the solar image. The standard angular resolution used 
in the regular observing programs (magnetograms of the large-scale magnetic 
fields, LSMF) is 100" and the scanning step is 91". To cover both lines with 
a single exposure, we used the fourth order of the spectrograph, where more 
than 2.7 nm are available. Part of the observations obtained with this dispersion 
was done with the standard spatial resolution ( Figures 7 and 8), but most 
with focused solar image and standard scanning step, what corresponds to a 
resolution of ~ 10". 

The cycle time (one phase of the electro-optical modulator - KD*P crystal 
in our case) was 30 ms, the signal integration time in every point of the scanning 
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process across the Sun was four seconds with and four seconds without the half- 
wave plate. The A/2 plate is periodically inserted into the light beam in front of 
the coelostat to monitor the zero-level position. 

An example of Stokes / and Stokes V/Ic profiles across the whole CCD 
detector is shown in Figure 1. The data correspond to the point next to the 
solar disk center, observed on 3 February, 2009 with a resolution of about 10". 
It is easy to see how much information is provided by such observations. 

To derive quantitative information about magnetic and thermodynamic pa- 
rameters from Stokes profiles such as those shown in Figure 1 is a complicated 
task. We use the approach (see Demidov et ai, 2002), which numerically imitates 
the work of solar magnetographs in the center-of-gravity mode (see Solanki, 
1993), to obtain the magnetic flux density (in the following we will use the term 
"magnetic field strength" for it throughout the paper) from the recorded data. 
The code calculates the position of the spectral line in two states (±A/4) of the 
polarization analyzer. Parameters for the "exit slits" are selected corresponding 
to the different lines. This position corresponds to the middle between the "slits" , 
where the intensities are equal. If Ai is the wavelength of the spectral line in one 
analyser state and A2 in the other one, the value AA = (Ai — A2)/2 is considered 
as the Zeeman shift associated with the magnetic-field strength B by the known 
equation AA = 4.67 x lO^^^A^i?, where B is in Gauss and A in cm. For the 
case shown in Figure 1 , the magnetic strength calculated this way is -26 G for 
Fei 525.0 nm and -64 G for Fei 523.3 nm. Parameters for the "exit slits" are: 
the width W is equal to 9.32 pm, the separation A (distance between line center 
and the "slits" centers) is equal to 6.21 pm for the first line and W — 24.84 pm, 
A = 15.53 pm for the second line. 

The data used for Figure 1 offer the possibility to estimate the errors of 
our measurements. To do this the noise level of the Stokes V/Ic signal in the 
continuum can be used. The best choice for it is the part of the spectrum on the 
right of Fei 524.4 nm. Calculations made for the ?» 150pm wide band show that 
the rms value in this region is 1.2 x 10~*. Taking into account that V/Ic =0.01 
corresponds to ~ 26 G for Fe i 525.0 nm and « 64 G for Fe i 523.3 nm, we conclude 
that the formal error of the measurements is w 0.3 G for the first line and « 0.8 G 
for the second. 

Figure 2 shows the Stokes / and Stokes V/ Ic profiles for this pair of lines from 
the same observations as in Figure 1 . Figure 3 compares the Stokes V/Ic profiles 
of both lines. Despite of the large difference in the (7-factors, the amplitudes 
of the Stokes y//c-profiles are almost the same and the peak positions of the 
Stokes V/Ic profiles are close to the steepest parts of the corresponding Stokes 
/ profiles. As a consequence of the fact that Fei 523.3 nm is deeper and wider 
than Fei 525.0 nm, the separation of the Stokes V/Ic extrema for the first line 
is almost two times larger than for the second one, although the splitting factor 
of the latter is larger. Therefore, in the case of weak fields, the separation of the 
Stokes V/Ic profiles peaks cannot be a measure of the magnetic field strength 
when lines with different widths are used. Indeed, if the magnetic-field strength 
is far from the Zeeman saturation regime, the locations of the Stokes V/Ic profile 
peaks are insensitive to the field and determined from the derivative of the Stokes 
/ profile. 



spl015.tex; 19/08/2009; 12:25; p. 4 



Solar Magnetic fields and the SOHO/MDI Calibration 



1.0 



03 
+^ 

13 

(D 
.> 

V- ' 

o 0.5 



CO 
(D 

O 

CO 





IT 



0.0 
522.8 



523.7 524.6 
Wavelength [nm] 



525.5 





0.010 




0.005 


o 




> 


0.000 


CO 




® 




O 


-0.005 


tn 






-0.010 




-0.015 




522.8 



523.7 524.6 
Wavelength [nm] 



525.5 



Figure 1. Spectra of Stokes parameters J (top panel) and Vjlc (bottom panel) along the 
whole wavelength band including Fel 523.3 nm and Pel 525.0 nm, covered by the detector. 
These observations were obtained close to the solar disk center with a magnetic-field strength 
of -26 G measured in Fel 525.0 nm. 
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Figure 2. Examples of Stokes / (in relative units and normalized to fit the scale of graph) 
and Stokes Yjlc profiles for Fel 523.3 nm (top panel) and for Fel 525. Onm (bottom panel). 
The same observations as in Figure 1 were used. 
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Already Figure 3 allows for the disk center a preliminary estimate of the ratio 
R = 5(523. 3)/B(525.0). In agreement with Demidov et al. (2008), we prefer the 
notation R for the ratio or the regression coefficient. The ratio of the Stokes V/Ic 
amplitudes (average of the amplitudes of the blue and red Stokes V/ Ic peaks) of 
the two lines is l/(523.3)/y(525.0) = 0.0097/0.0106 = 0.915. Normalizing to the 
corresponding Lande factors, we obtain R = 0.915x (3.0/1.3) = 0.915x2.3 = 2.1. 
This value is rather close to the estimation R = B(523.3)/B(525.0) = (-64 G/- 
26 G) = 2.4, obtained from calculations of the magnetic-field strengths using 
wide parts of the line wings. 

One point is not enough for a reliable estimation of the 5(523. 3)/i3(525.0) 
ratio. The scatter plot, calculated for three full-disk magnetograms observed on 
1, 2, and 3 February, 2009, is shown in Figure 4. Parameters for the "exit slits" 
are the same as mentioned above. We see a rather high correlation coefficient 
p = 0.93, and a systematic difference with regression coefficient R = 1.97 ±0.02. 
To calculate i?, the reduced major axis method (Davis, 1986) is used. A formula 
to estimate the error (which is determined from the scattering of the points, and 
which is less when the correlation coefficient is higher) is taken from this book 
too. According to our measurements, the average factor that should be used to 
adjust measurements in Fei 525.0 nm to those in Fei 523.3 nm, is 1.97 ± 0.02. 

The next logical step is to check for possible spatial variations of the ratio 
R = -8(523. 3)/i?(525.0) across the solar disk. For this purpose, we divide the 
solar disk into polar and equatorial sectors as in Demidov et al. (2008). We derive 
a polynomial fit R = 1.74 — 2.43^-1-3.43//^, shown in Figure 5. There are no 
significant differences in the CLV of R between polar and equatorial sectors of 
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Figure 4. Correlation and regression analysis of solar magnetic field observations in 
Fe I 523.3 nm and Fe I 525.0 nm. Data for three days of observations, 1, 2 and 3 February, 2009, 
are used. is the number of points, p is the correlation coefficient, and R is the coefficient of 
linear regression (the slope of the line through the scatter plot). 

the disk. At disk center, R is 2.74, and closer to the hmb, at /i = 0.3, R is 1.32. 
The dashed hne is a quadratic polynomial fit through the observed points. 

A significant decrease of the ratio i3(523.3)/i?(525.0) with increasing helio- 
centric distance is illustrated in Figure 6, where V/I^ profiles of the analyzed 
lines are shown for the point next to the south pole with fi = 0.56 and with 
magnetic field strengths B(525.0) = 11.6 G and B(523.3) = 19.8 G. The ratio of 
these strengths, R ~ 1.5, is rather close to i? = 1.45 from the analytical formula 
in Figure 5 at ^ = 0.56. 

As mentioned above, observations in Fei 523.3 nm and Fei 525.0 nm with 
the STOP traditional spatial resolution of 100" were performed besides the 
observations with the rather high spatial resolution in order to look for a possible 
dependence oi R = B(523.3)/S(525.0) on the spatial resolution. Ulrich (1992) 
found a significant increase of R with decreasing spatial resolution: R = 4.5 at 
disk center at a resolution of 20" instead of i? = 3.9 at a resolution of 5". 

Examples of full-disk magnetograms with a resolution of 100", obtained on 1 
February, 2009 are shown in Figure 7. Despite of the very weak magnetic field 
strengths everywhere across the disk (minimum of solar activity, there were no 
sunspots and active regions), both maps show identical spatial structures, and 
they are well correlated. This is proved by the regression and correlation analysis 
presented in Figure 8. Calculations of the CLV of i? = B(523.3)/B(525.0) for 
these data yield: R = 2.24 0.03^ -^0.56^^. We see that R is 2.83 at disk center, 
slightly more than the value of i? = 2.74 obtained with high-spatial resolution. 
Despite the decrease of the resolution by a factor of ten. we encounter only a 
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Figure 5. Center-to-limb variation of the ratio R = S(523.3)/B(525.0) for the polar (N-S) 
and equatorial (E- W) sectors of the solar disk. 
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Figure 6. Example of Stokes Yjlc profiles of Fel 523.3 nm and Fel 525.0 nm for a point near 
the south pole with ^l = 0.56 and with a magnetic-field strength B = 11.6 G from Fel 525.0 nm. 
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Figure 7. Full-disk magnctograms observed with a spatial resolution of 100" in Fe I 523.3 nm 
(left panel) and Fcl 525.0 nm (right panel) on 1 February, 2009. 




-4 -2 2 4 

B525.0 [G] 

Figure 8. Correlation and regression analysis of solar magnetic field observations in 
Fel 523.3 nm and Fcl 525.0 nm for observations with a spatial resolution of 100", made on 
1 February, 2009. Notations arc the same as in Figure 4. 

small « 3% increase of R. The difference is much larger for large heliocentric 
distances, but most probably it is caused by a weakening of the field strengths 
and by an increasing influence of the noise. 

3. Discussion and Conclusion 

The cornerstone of the study by Ulrich et al. (2009) is the difference of magnetic 
field strengths measured in the different parts of the profile of Fci 523.3 nm, i.e. 
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the sampling of the hne. They wrote: "contrary to the assumptions of Frazier 
and Stenflo (1972), Stenfio (1973) and Frazier and Stenflo (1978) we find that the 
field indicated by A 5233 A depends on the spectral sampling" and that "we [...] 
recommend the use of a sampling point as close to the line core as is practical" . 
In the work of Ulrich (1992). the parameters of the spectral configurations with 
fiber-optic image reformattors were: centered on Fei 525.0 nm - separation A = 
±3.9 pm; centered on Fei 523.3 nm: - A = ±4.5 pm. 

In the recent work, Ulrich et al. (2009) obtained new MWO observations 
with the same sampling for Fei 525.0 nm and with a new set of separations of 
±0.9, ±2.9, ±8.4, ±10.2,and ±17.7pm for Fei 523.3 nm. Computations of mag- 
netic fields from various spectral sampling pairs and at ±8.4 pm have shown 
different values of correlations coefficients p and slopes in the corresponding 
scatter plots, depending on the value of the separation and the range of 9 (see 
Figures 1 to 3 and Table I in Ulrich et al., 2009). The correlation coefficients 
for measurements at ±8.4 pm with the other ones were: p = 0.50 at ±0.9 pm, 
p = 0.90 at ±2.9 pm, p = 0.90 at ±10.2 pm and p = 0.50 at ±17.7 pm. Only 
the correlation coefficients for the pairs nearest to ±8.4 pm are reliable. In two 
other cases they arc too small for a solid statistical analysis. The examples of 
slopes obtained by Ulrich et al. (2009) for the central area of the disk (range 
of sine from 0.0 to 0.4) are: 1.665 ±0.067 at ±0.9 pm, 1.795 ±0.027 at ±2.9 pm, 
1.162 ±0.010 at ±10.2 pm and 2.694 ±0.094 at ±17.7 pm. 

The fact that magnetic fields measured in points near the line center (±2.9 pm)| 
are weaker (by a factor of 0.73 after application of additional assumptions and 
calculations, see Figure 12 of Ulrich et al., 2009) than those measured in a point 
near the middle of the line wings (±8.4 pm) is the basic argument for the decrease 
of the correction coefficient i?, following from direct comparison of data sets in 
Fei 523.3 nm (sampling ±8.4 pm) and Fei 525.0 nm. For the disk center, the 
value R = 4.12 is obtained instead of 5.5. 

We are convinced that only the Stokes V/ Ic profile over the whole spectral line 
can characterize the magnetic field properties in the range of the line formation. 
But for the comparison with MWO results, it is important to find relationships 
between strengths in the different parts of Fei 523.3 nm in our data. For this 
purpose, calculations with different "exit slit" parameters, close to those recently 
used at MWO, were performed. The results for 3 February, 2009 with 10" spatial 
resolution are given in the following. 

The values of "slit" widths and separations are expressed in the integer num- 
bers of CCD pixels. In observations, the width of one pixel is equal to 0.777 pm. 
Following Ulrich et al. (2009), all data were analyzed relatively to the measure- 
ments in the part of the wings with a separation of A = 8.5 pm and with a "slit" 
width oi W = 9.32 pm. The comparison of data with these "slit" parameters 
with other ones are given in Table 1. For the points at ±3.1 pm, our result 
differs from that of Ulrich et al. (2009). According to our data, magnetic fields 
measured closer to the line center are almost the same (only 1.07 times weaker) 
as measured at ±8.5 pm. That means, if we correct measurements at ±8.5 pm 
taking into account the regression coefficient i?(±8.5 pm)/i3(±3.1 pm) = 1.07, 
we get a small difference to measurements in Fei 525.0 nm. Indeed, Ulrich et 
al. (2009) found 5(523.3 ± 8.4pm)/B(525.0 ± 3.9 pm) = 5.5 for the disk center. 
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Table 1. Comparison with difTcrent "slit" parameters. 



AA [pm] Slit width W [pm] p _R, = B(±8.5pm)/B(±AA) 



±15.53 


24.82 


0.991 


1.271 ± 0.009 


±0.8 


0.80 


0.870 


1.018 ± 0.028 


±3.1 


1.55 


0.977 


1.073 ± 0.013 


±10.1 


9.32 


0.998 


1.027 ± 0.003 


±17.8 


9.32 


0.935 


1.528 ± 0.030 



Applying a special correction factor, they obtained the value of 4.12. In our case 
we find 5.5 x (1/1.07) = 5.1. 

All results described above have a direct impact on the problem of the calibra- 
tion of SOHO/MDI magnetic- field data. Some correction coefficients suggested 
by MWO scientists were already mentioned before. The correction factors for 
the MDI data are much less if we use the values derived from our observations 
instead of those obtained by Tran et al. (2005) and by Ulrich et al. (2009). For 
disk center, we obtain 1.61 /(4. 15/2.74) ~ 1.06 instead of 1.61 according to 
Ulrich et al. (2009). 

Our independent estimations of the SOHO/MDI correction factors are ob- 
tained using our previous results from the comparison of SOHO /MDI data with 
SSO observations (important: both analyzed data sets correspond to spatial 
resolution 100"), published by Demidov et al. (2008). From Figure 1(e) of that 
paper we get an average value of 2.75 for the B(SOHO/MDI)/B(SSO) ratio. In 
order to correct the SOHO/MDI data (when they coincide with the corrected 
data in Fei 525. Onm), we have to multiply them with the coefficient 1.97/2.75 
= 0.72, if we take R = 5(523. 3)/S(525.0) = 1.97 from our data with 10" spatial 
resolution, and 2.55/2.75 = 0.93, if we take R = 2.55 from data (see Figure 8) 
with 100" spatial resolution. Our correction coefficient is smaller by a factor of 
about two compared to the recent MWO results (Tran et al., 2005, Ulrich et al., 
2009), and instead of increasing the magnetic field strengths from SOHO/MDI 
data we have to decrease them. We would like to note here, that, according to 
the study of Demidov et al. (2008), observations at Mount Wilson and Sayan 
observatories made in Fei 525.0 nm are in almost perfect agreement with each 
other. Therefore, the problem lies in the differences between the observations in 
the strong line Fei 523.3 nm. 

It seems that this result contradicts to the result of Berger and Lites (2003), 
who obtained a correction factor of 1.56, based on a comparison of SOHO/MDI 
and Advanced Stokes Polarimeter (ASP) data. But observations with the ASP 
were made in a different spectral line (Fei 630.25 nm), and using different lines 
can lead to quite different results (see Gopasyuk et al., 1972 and Demidov et 
al., 2008). Further, Berger and Lites (2003) analyzed only observations made in 
active regions, instead of full-disk quiet-Sun observations in our case. Besides 
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that, the observations with ASP and SOHO/MDI were made with much higher 
spatial resohition. 

In connection with this issue it is important to note that comparisons of 
SOHO/MDI with Kitt Peak spectro-polarimeter data (Jones and Ceja, 2001, 
Wenzlcr, Solanki, and Krivova, 2005) caused a decrease of the SOHO/MDI 
measurements by a factor of « 1.4 to adjust them to the Kitt Peak data. 

Summarizing, we arrive at the conclusion that the issue of calibrating SOHO/ 
MDI data is rather far from being solved, and new investigations and observa- 
tions are urgently needed. 
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